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Abstract 

The persistent current in three-dimensional mesoscopic rings is investigated 
numerically. The model is tight-binding one with random site-energies and 
interaction between electrons. The Hartree-Fock approximation is adopted 
for the interaction between electrons, and models with up to 3 x 4 x 40 sites 
are investigated. It is shown that the long-range Coulomb interaction en- 
hances the persistent current considerably for rings with finite width. It is 
also shown that long-rangedness of the interaction is essential for the enhance- 
ment. Screening of the random site-energy is attributed to this enhancement. 
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It has been predicted by Biittiker et al. [T[| that an equilibrium persistent current flows 
in a mesoscopic normal- metal ring threaded by a magnetic flux. This prediction has been 
experimentally confirmed 10-0. Especially, the experiment on single Au ring || revealed 
that the size of the persistent current observed was quite large. Namely the observed current 
was I = (0.3 ~ 2.0)/o, where Jo = evp/L is of the order of the persistent current in an ideal 
one-dimensional ring with vf being the Fermi velocity and L being the circumference of the 
ring. The current is periodic in the magnetic flux piercing the ring with period 0o — h/e 
as expected. Thus the persistent current was predicted theoretically and confirmed by the 
experiment. However, this is not the whole story. Quantitative understanding of the current 
has not been established. 

It has been clarified by numerical 0|| and theoretical investigations that theory 
without taking into account the interaction between electrons can give current which is too 
small compared to the experiments 0,[§. Inclusion of the mutual interaction into theory is 
hard and difficult to assess the validity of the approximation employed, but several groups 
|||5|] have done such calculations, and obtained results that sample averaged current is of the 
order of (l e /L)I . In this case the current is periodic in 4>o/2, because sample averaging kills 
the component with period 0o- This current is much smaller than J in the experimental 
situation of L » l e , where l e is the elastic mean free path. In these theories the Coulomb 
interaction acts mainly to maintain the local charge neutrality of the system. Another 
theory by Kopietz takes into account long-range Coulomb interaction energy associated 
with the charge fluctuation. This theory gives the average persistent current of the order of 
O.l/o with period </>o/2. This value is large enough to explain the experiments |2j. However, 
since the fluctuation should be suppressed by the Coulomb interaction itself to maintain 
the approximate local charge neutrality, it is not yet certain that such large value survives 
improvement of the theory. 

On the other hand, numerical investigation is free from approximations, and has a pos- 
sibility to give a clear answer to the origin of the large persistent current. However, unfor- 
tunately rigorous treatments are limited to quite small system sizes. Thus until now only 
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small one- dimensional systems have been investigated with negative answer to the role of 
the Coulomb interaction |p"T|-|i"3|j. Namely, it has been clarified that in one- dimensional sys- 



tems the repulsive interaction between electrons suppress the persistent current in situations 
relevant to the experiment. We have also investigated the one-dimensional system with the 



Coulomb interaction between electrons \TA \ . Our results were the same as others. However, 
we also found that the Hartree-Fock approximation gives almost the same results as those 
of the exact diagonalization. The merit of the Hartree-Fock approximation is that we can 
deal with much larger systems. Thus we now apply our method to three-dimensional rings 
to investigate the effect of the Coulomb interaction. 

In the present work we have found that for the three-dimensional ring the Coulomb in- 
teraction enhances the persistent current. We have clarified conditions for the enhancement. 
We have found that long-rangedness of the Coulomb interaction is essential to this enhance- 
ment: Short-ranged interaction does not enhance the current. We also found a condition for 
the cross section of the ring to have the enhancement. 

We consider the Anderson model with Coulomb interaction between electrons. For sim- 
plicity we neglect spin freedom of electrons. Thus we consider a tight binding model on a 
cubic lattice. Our model ring consists of N# x N w x sites: circumference N^a, width N w a, 
and height Nha, where a is the lattice constant. The coordinate of the z-th site is given by 
three integers (£i, Wi, hi), where 1 < li < Ng, 1 < Wi < N w , and 1 < hi < N^. Each site has 
random site-energy which has uniform distribution between — W/2 and W/2. Thus our 
Hamiltonian is 

H = - y,(ti.j<'" '< : i<-j + h.c.) + 5> iC jci 

+ \ V —( c Ui - Po)(c]cj - po), (1) 

where cj(cj) is the creation (destruction) operator of a spinless electron, tij is non-zero and 
equal to t only for nearest neighbor hopping, dij is non-zero only for hopping along the 
circumference of the ring and has a value 2ir<f>/Nt<f)o, V = e 2 /47rea is the strength of the 
Coulomb interaction, and po is neutralizing positive charge at each site. To make a ring 



which has the same lattice constant everywhere we have embedded our three-dimensional 
ring in the four- dimensional space, so the distance between i, j sites is given by 

(^) 2 = 2i? 2 {l-cos[^ -*,)]} 

+ ( Wl - w 3 ) 2 + (hi - hj) 2 , (2) 

where R = [2 sin(7r/i\^)] _1 is the radius of the ring. 



The following treatment is the same as before fL4| : The interaction term is decoupled 
following standard Hartree-Fock scheme by introducing the average (c\cj), which is deter- 
mined self-consistently by iteration. Once the ground state energy E g ((f>) is determined 
self-consistently, the persistent current 1 (0) is given by the derivative of the energy E g ((j)) 
with respect to the magnetic flux: 

Since E g (<f)) is an even function of and periodic with period 0o, we first Fourier analyze 
the energy, 

E{(j>) = E + J2E n cos(27m^-) , (4) 

and obtain the Fourier components of the current: 

T . , . J2, 27m „ . . 

n=l 00 00 
°° t 

= y2l n —sm(27m—). (5) 

£i 0o 0o 

Here is a dimensionless number, which depends on realization of the random site-energy 
and the strength of the Coulomb interaction. In the following we mainly concentrate on the 
0o-periodic component, I\. 

We have investigated up to 480-site systems by this method, the largest being Nt = 
40, N w = 4, and Nh = 3 with 200 electrons. The randomness is chosen such that the 
sample is in the diffusive regime, which is realized around W/t = 2 for the samples we 
have considered. The results for the largest four samples are shown in Figs.l, 2 and 3. 



First in Fig.l dependence on the randomness is shown for the non-interacting case to clarify 
the situation before we take the interaction into account. An arrow indicates the value 
of I\ = 4k /Ni. The size of this current (4ftfNg)t<j)Q is equal to cvf/L = I , if we use 
vf = 2\t\a/h; namely it indicates the size of the persistent current for a pure one-dimensional 
ring of the same circumference. Figure 2 shows how the interaction enhances the current. It 
is seen that the interaction enhances the current beyond the value of the pure 1-d system, 
Jo- On the other hand Fig.3 shows the results obtained by replacing the long-range Coulomb 
interaction by short ranged one, namely V x (a/r^-) is replaced by VS{rij — a). It has the 
same value for the nearest neighbor but vanishes otherwise. The results show essentially no 
enhancement for this case. This comparison clarifies the origin of the enhancement: long- 
rangedness is essential for the enhancement. Quite similar results have been obtained for 
smaller samples (33 x 3 x 2 sites and 16 x 3 x 2 sites). It should be noticed that the sign of 
l\ is not fixed: In these figures the positive values are shown by open symbols and negative 
ones are shown by closed symbols. It can be both positive (paramagnetic) or negative 
(diamagnetic) depending on samples. 

Since our present result is opposite to that of one- dimensional system, we investigated 
the boundary between the enhancement and the suppression. Both two-dimensional and 
tree-dimensional rings were examined. For the two-dimensional system, N w was kept to 
unity and was increased from one to six. Definite enhancement by V/t was observed 
for systems with > 4. For the three-dimensional system, N w was kept to two and is 
increased from one to three. It was found that > 2 is enough to have the enhancement. 
Thus in both cases enhancement occurs for systems with N w x > 4; otherwise the current 
is suppressed. 

In the present work we have shown for the first time that the Coulomb interaction can 
enhance the persistent current beyond the value for an ideal one-dimensional ring, I |L5] . 



Then the next question will be if the parameters used here are relevant to the experimental 
situation or not. We think that choice of the randomness W/t poses no question. We need 
to be in the diffusive regime, and within this regime V/t dependence is qualitatively similar. 



Appropriate choice of V/t is a difficult question, since size of our system is still too small 
compared to the experiment, where about 10 8 atoms are involved. Here we only comment 
that if our system is actually made of Au atoms, t is of the order of leV, and bare Coulomb 
energy between one lattice constant is larger than t: V — e 2 /47rea ~ 5eV > t. So even if 
we consider the optical dielectric constant e(oo) by ion core, V and t are of the same order. 
We have found strong enhancement even when V is one order of magnitude smaller than 
t. In the actual system the interaction is screened by electrons, but such screening effect is 
included in our self-consistent solution. We need the flux dependence of the total ground 
state energy which includes the energy associated with the screening. Therefore use of the 
long-range Coulomb interaction is needed and justified. 

Our treatment is not exact as we use the Hartree-Fock approximation. We justify it in 



two ways. Firstly, as shown in ref. [[L4[] the results by our Hartree-Fock approximation is con- 
sistent with those by the exact diagonalization for small one-dimensional systems. Secondly, 
in this work we have clarified that the long-range part of the interaction is essential for the 
enhancement. It is expected that correlation effects beyond the Hartree-Fock approximation 
mainly affects the short-range behavior. However, since we consider spinless electrons, main 
part of the short-range correlation, on-site correlation, has already been taken into account. 
Thus we expect that our present result is almost the same as what would be obtained if the 
correlation were taken into account exactly. 

The final question will be the mechanism of the enhancement. One possible expla- 
nation is the screening of the random potential. The electron at the Fermi level feels 
screened site-energy which is flattened, so the current is enhanced. In the case of one- 
dimensional system the screening is also effective, but electrons are not allowed to pass 
through each other, thus the current is suppressed by localization of some of electrons. 
When the cross section of the ring is small enough (e.g. N w x < 4), this block- 
ade is still effective. This explanation is justified as follows. We have examined which 
term in the Hamiltonian mainly contribute to the current. We found that the contribu- 
tion from the first two terms of the Hamiltonian Eq. (|l|), i.e. the single-particle terms, 
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I s (<j>) = —(d/d4>)(—Y,i^j(ti,je l9ij clcj + h.c.) + Y^i^icjci), amounts to about 90% of the to- 
tal current, where ( ) means expectation values with respect to the self-consistent ground 
state. The remaining 10% of the current comes from the interaction term of the Hamiltonian 
which consists of the Hartree contribution, Ih{4>) = —(d/d(j)) Yltyj V{ a /2 r ij){ c l c i){ c ] c j) , and 
the Fock contribution, = (d/d<j)) J2i^j ^ / (°/2^ij)(cJcj)(cjcj), the latter being usually 

smaller than the former. Thus the enhancement is caused mainly by the single particle part 
of the Hamiltonian. This means that the ground state is quite different from that of the 
non- interacting system, where we get quite small persistent current. We expect that the 
electron states near the Fermi level are more like those of the system with W = V = 0. 
This explains why the theories using the diagrammatic expansion P,p, l6|,17| could not get 
the enhancement: Since the ground state has changed considerably from that without the 
Coulomb interaction, perturbational calculation cannot give correct answer. 

In conclusion we have shown for the first time that the persistent current is enhanced 
beyond the ideal one- dimensional value by the long-range Coulomb interaction for three- 
dimensional rings. We have clarified that a ring with finite cross section and long-rangedness 
of the interaction are necessary for the enhancement. We attribute the enhancement to the 
screening of the random site-energies. 

One of the authors (D. Y.) thanks Aspen Center for Physics where part of this work was 
done. The numerical calculation was done by HITAC S-3800/480 at the Computer Center 
of the University of Tokyo. 
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FIGURES 

FIG. 1. Fourier component of the persistent current with period 4>q is plotted as a function of 
the site randomness W/t for the non-interacting system (V/t = 0). The four different symbols, 
circle, two kinds of triangles, and square, indicate four different realization of the randomness. The 
absolute values of the current are plotted with the open (closed) symbols indicating that the value 
is positive (negative). An arrow with letter Iq indicates the value of I\ where Iit/fo = Iq. 

FIG. 2. Fourier component of the persistent current with period 4>q is plotted as a function of 
the strength of the Coulomb interaction V/t for the randomness W/t = 2.5. The same four samples 
in Fig.l are used with the same symbols. Here also the open (closed) symbols mean that the value 
is positive (negative). An arrow with letter Iq indicates the value of I\ where ht/fo = I®. 

FIG. 3. Fourier component of the persistent current with period 4>o is plotted as a function 
of the strength of the short-range interaction V/t for the randomness W/t = 2.5. The same four 
samples in Figs.l and 2 are used with the same symbols. Here also the open (closed) symbols mean 
that the value is positive (negative). 
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